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Interaction energy surfaces for the interaction of Li(2 2 S) and Li(2'-P) with H 2 are calculated using 
approximate Hartree-Fock trial functions. The crossing of energy curves is observed for C2 V confor- 
mations if the H 2 internuclear distance is sufficiently large. No crossing is observed for colinear col- 
lissions for any H 2 distance. 

The wave function of the strongly attractive state involved in the crossing is related to the metasta- 
ble negative-ion states that are postulated to account for resonant electron-molecule scattering. Such 
a charge-transfer state can only be bound for C 2v conformations for the H 2 molecule. The likelihood 
and geometry of the crossing complexes for other molecules interacting with alkalis is discussed in 
terms of the formation of these resonance charge-transfer states. 
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1. Introduction 

Fundamental to the understanding of electronic- 
vibrational energy transfer involving alkali atoms 
and diatomic molecules is an accurate representa- 
tion of the energy of interaction. The qualitative 
features of these surfaces was first considered with 
semi-empirical techniques [l]. 1 To this date there 
are no ab initio calculations. Accurate energy cal- 
culations utilizing correlated wave functions are 
not yet possible but these systems can be investi- 
gated immediately with available Hartree-Fock 
(H.F.) programs [2]. This is of value for two reasons. 
First, the approximate H.F. energy surface is con- 
siderably more reliable than the earlier results to 
the extent that the qualitative features of the sur- 
faces can be fixed. In particular, the question of the 
angular distribution of the interaction energy can 
be settled and an estimate given given of the acti- 
vation energy, if any, of the excitation of alkali D 
radiation by vibrationally excited diatomic s. Second, 
the molecular orbital description is more natural 
to the discussion of these reactions in the light of 
recent studies of electron-diatomic molecule resonant 
scattering [3]. 

This note will examine the qualitative behavior 
of the Li(2 2 S) and Li(2 2 P) interaction with the H 2 
molecule without chemical reaction as the simplest 
example of the alkali-diatomic system. There is a 
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total of four curves, two repulsive ones, and one 
weakly, and one strongly attractive curve which are 
depicted in figure 1. Except for the strongly attractive 
state there is little overlap of the alkali and diatomic 
charge distributions. It is our contention that the H.F. 
surfaces are qualitatively correct and of semiquanti- 
tative accuracy for this type of surface. The near con- 
stancy of the correlation energy over a wide range 
of geometries is presumed here [4]. 

The correlation situation is quite different for 
the strongly attractive curve. This curve has long 
been characterized as a charge transfer state [1]. 
The H.F. function does, in fact, describe a charge 
transfer state and is related to the ionic valence 
bond configurations. The negative ion H 2 state 
formed here is also related to the lowest energy 
resonance state of H 2 considered responsible for 
enhanced vibrational excitation at low incident 
electron energies [3]. Three electrons are now in- 
terpenetrating and the correlation [5] is similar 
to that in HeJ and H 3 . However, there are offsetting 
H.F. and correlation errors that are strongly depend- 
ent on the geometry. The depth of the H.F. attractive 
well should approximate the correlated one. 

Since even an accurate H.F. calculation of the 
strongly attractive curve can only be considered 
qualitative, we have contented ourselves with self- 
consistent-field solutions far from the H.F. limit. 
More accurate solutions will be considered when they 
are needed as a base for a correlated calculation. 
The H.F. techniques are briefly described and the 
results for a rather small but meaningful basis are 
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Figure la. Interaction energies for Li( 2 S and 2 P) with the ground state ofH 2 for C 2v 
geometries as a function of the H 2 internuclear distance. 
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FIGURE lb. Interaction energies for Li( 2 S and 2 P) with the ground 
state of rl 2 for C xv geometries as a function of the H 2 internuclear 
distance 

Energies in a.u. 



given in the next section. Similar calculations were 
completed for the isoelectronic LiHe system at the 
same level of accuracy. This permits some calibra- 



tion of the approximate H.F. values since limited 
configuration interaction results [6] exist for LiHe. 
In the final section the H.F. results are compared 
to the earlier semiempirical calculations. The angu- 
lar distribution and the position of the crossing in 
the interaction curves are particularly examined and 
earlier deductions on the linear surface are ques- 
tioned. Comparison of the H.F. wave functions for 
the charge-transfer state to the assumed resonance 
states in electron scattering is pursued to the point 
of speculation with regard to other systems. 

2. Hartree-Fock Calculation 

For doublet systems with one open-shell orbital 
outside of a closed shell core, the H.F. solutions 
satisfy a variant of the Brillouin theorem that was 
noted by Davidson [7]. If the vectors are designated 
Ci, a where ex designates the open-shell vector for 
which self-consistency is achieved then 



and 



v>i, a^^j, a bij 
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where S is the basis function overlap matrix and H 
is the total electronic Hamiltonian. ^P*, a is the H.F. 
function where the excited (fi orbital replaces the 
<p a orbital in M 7 ,*, a . (Calculation of any one state will 
yield upper bounds to higher Rydberg states. The 
accuracy of the bounds depends upon the size of the 
basis set and its ability to fit the higher Rydbergs and 
the extent to which the core orbitals relax for differ- 
ent ^ tt , a - As noted in the introduction no special 
effort will be made to achieve the H.F. limit. But 
a reasonably large set must be used to insure correct 
qualitative behavior. Both the Li( 2 S) and H( 2 S) bases 
are chosen from the smallest of the accurate Gaussian- 
type functions (GTF) sets given by Huzinaga [8]. 
The 2s H , 2p Li , 2p H , 3s Li , and 3p Li orbitals are approxi- 
mated by choosing effective Slater exponents for the 
appropriate atomic Rydberg orbitals and using the 
large coefficient Gaussian fits to the Slater orbitals [9], 
This procedure is very rough since not all the GTF 
required can be used because of computer time lim- 
itations so that compromise exponents must be chosen. 
This set has been contracted severely with the re- 
sultant H 2 and Li energies in table 1. The Li( 2 S- 2 P) 
interval is 0.11 eV larger than the H.F. value and 0.12 
eV larger than the experimental value [10]. The 
error is almost entirely in the 2s Li due to the use of 
ls Li coefficients to severely contract the s u set. For 
H 2 the error in the H.F. energy is about 0.6 eV at 
1.5 a.u. compared to an accurate H.F. energy. This 
large discrepancy is due to the contraction of the 
H 2 basis using the Huzinaga expansion coefficients 
without scaling the orbital exponents. Since the con- 
tracted function is equivalent to a single Slater- 
type function (STF) the Is function exponents de- 
termined by Fraga and Ransil [11] can be used as a 
guide. The larger exponent would improve the en- 
ergy significantly but would effect improvements 
mainly in the H 2 binding and not in the interaction 
region between the Li and H 2 which is the main 
concern here. More flexible bases will be used when 
the correlation is calculated in subsequent work. 
The errors due to the small basis considerably ex- 
ceed those that arise from the use of orbital energies 
to compute the total energies. The general shapes 



Table 1. Energies of asymptotic atoms and molecule 



of all curves can be found in figure 1 for the neutral 
states and in figure 2 for the Li + + H 2 interaction. 





-E 


• 


Atom/molecule 


Present calculation 


Accurate H.F. 


He 


2.86113 


"2.86168 


Li + 


7.23597 


11 7.23641 


Li(*S) 


7.42714 


c 7.432724 


Li(*P) 


7.36342 


c 7.365068 


H 2 (1.0) d 


1 .02984 




H 2 (1.5) 


1.11010 




H,(2.0) 


1.08197 




H 2 (2.5) 


1.03299 




H,(1.4) 




* 1.13354 



"Energy in atomic units: a.u. of energy =e*/oo= 27.210 eV. 
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Figure 2. Interaction energies for Li^S) with the ground state 
of H 2 for (a) C% v and (b) C xV geometries as a function of the H 2 
internuclear distance. 



energies in a.u. 



The quantitative value of the B> curves is limited 
by an uncertain correlation error. Correlation effects 
will be considered in two categories. The first type 
are terms which correlate the asymptotic fragments. 
The K shell of the Li atom is, of course, essentially 
unaffected by the interaction. It is presumed that 
for all but the B 2 state the correlation of the H 2 bond 
is only slightly changed for moderate interaction 
energies. The correlation energy curve for H 2 is quite 
constant to about 1.5 a.u. and then starts to rise 
rapidly [14]. For the B 2 state the separation of the two 
hydrogen atoms in the presence of the Li atom does 
not lead to a large increase in the correlation energy 
as /?(H — H) exceeds 1.5 a.u. Since a 6 2 (lo-„) orbital 
is occupied, the asymptotic fragments would be 
H+H with Li for electrical neutrality. The cor- 
relation energy [15] for H is only a little less than 
that for H 2 and the usual large asymptotic correla- 
tion error is avoided. For H 2 separations less than 
about 1.5 a.u. the correlation energy in the H 2 pair 
should be comparable for the Ai, Bi, and B 2 states. 

At R(U — H) = 2.5 a.u. the B 2 state would be about 
0.01 a.u. lower in energy than the Ai and Bi states 
due to the difference in correlation energy in the H 2 
electron-pair for the two categories of states. Off- 
setting this error would be the increase in the sec- 
ond category of the correlation energy due to the 
three interpenetrating electrons in the B 2 state. For 
HeJ the correlation energy at equilibrium exceeds 
that for the He atom by about 0.025 a.u. It is not 
expected that the additional correlation energy in 
LiH 2 2 B 2 will be that large since the charge is not 
completely transferred from the Li atom in the range 
of conformations of interest. Assuming a direct pro- 
portionality of the excess correlation energy to the 
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amount of charge transferred, the value is 0.008 a.u. 
for the curve where 7?(H — H) is 1.5 a.u. On this basis 
it is felt that the calculated surfaces are semiquanti- 
tatively correct. For these surfaces the zero activa- 
tion energy crossing for activation of the 2 P state by 
colliding the 2 S atom with a vibrationally excited 
molecule occurs at a conformation where /?(H — H) 
is a little less than 2.0 a.u. 

An oblique check on tnese calculations can also 
be made by considering the isoelectronic diatomic 
interaction Li + He for which accurate calculations 
exist. The asymptotic error is still due to the poor 
Li(-S) function. Only the 2 n state was calculated di- 
rectly and the well depth and spectroscopic constants 
for the 2 I1 states are compared in table 2. All the 
LiHe interaction energy curves are exhibited in 
figure 3. 

Table 2. Spectroscopic constants and 
dissociation energy, LiHe 
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FIGURE 3. Interaction energies for Li( 2 S and 2 P) with the ground 
state of He. 
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3. Discussion 

The present results confirm the early Heitler- 
London predictions of crossings in the energy 
curves for triangular geometries [la]. There are 
no crossings in the colinear cases. In fact the linear 
curves strongly resemble the LiHe interaction en- 
ergy curve. Laidler [lb] has discussed the Na( 2 P) 
quenching in terms of the linear complex and has 



argued that the conclusions would not be affected by 
the shape of the complex. For the Li reactions the por- 
tion of the surface calculated here would argue 
quite the opposite. It is improbable that colinear col- 
lissions of vibrationally excited or dissociated H 2 
with Li( 2 S) would lead to excitation to Li( 2 P) or the 
converse. But crossings are available in the C 2v 
configuration and there is no classical activation 
energy. The curves indicate a very rapid change in 
crossing position and well depth as a function of the 
internuclear IT distance. This would permit the 
strong coupling of the electronic and vibrational de- 
grees of freedom which is in accord with both past [lb] 
and present [12] interpretations of these reactions. 
In particular deactivation of alkali resonance radia- 
tion will occur predominantly by vibrational excita- 
tion. 

Two features of the curves are worth specia^ 
mention. First, the crossing in the Hartree-Fock ap- 
proximation is a kind of conical intersection [13] of 
energy curves. The doublet open-shell configurations 
have all off-diagonal energy matrix elements zero 
even if the symmetries are identical. 2 A' states 
would exhibit a conical intersection in the H.F. ap- 
proximation. However, it is difficult to see how the 
degeneracy could be maintained when the full cor- 
relation interaction is introduced although the split- 
ting is probably small. The second feature is the al- 
most congruence of the H.F. 2 Bi or 2 n neutral curves 
and the ion curve for the same geometry. The well 
depth in the n-type states is not due, in this approxi- 
mation, to any slight binding in the it orbital but pre- 
dominately to the ion core interaction with a polar- 
izable, unshielded atom. Similar results are found 
for LiHe. Although the well depths are about 40 per- 
cent smaller than those of Schneiderman and Michels 
or an accurate H.F. result, the qualitative behavior is 
similar. The 2 n and ] X + well depths are quite com- 
parable with the 2 n dissociation energy greater by 
0.01 eV. 

The large H 2 separations required for the crossing 
were previously linked to the 2 2j state [1] of H 2 . 
H.F. wave functions can be similarly interpreted. For 
the B 2 state the p orbital is oriented along the H 2 
molecule. As the Li( 2 P) approaches H 2 , a fraction of 
the charge transfers to the H 2 in a charge distribu- 
tion very much like that of an isolated lcn,H 2 orbital. 
This is a penetrating valence orbital with strong 
antibonding characteristics which determines both the 
deep well and the large H 2 equilibrium separations. 
The picture that develops is that of a charge-transfer 
state arising from a low ionization potential atom and 
a diatomic molecule that can support a valencelike 
negative-ion resonance state. Those molecules which 
have no available low-lying valence-type orbitals 
will not have deep attractive interactions with excited 
alkali atoms. In this category would fall the Ne-iso- 
electronic molecules and saturated hydrocarbons. 
Quenching cross-sections are small for such mole- 
cules [16]. On the other hand the diatomics like 
N 2 have the iTy and cr u type valence orbitals avail- 
able. The n^ state can arise from a p orbital in a co- 
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linear collision for an orientation perpendicular to 
the line of centers. The C2 V geometry requires a d 
orbital on the alkali since the 7T g orbital correlates 
with an A 2 representation. But the a u antibonding 
valence-type orbital plays the same role in N 2 as it 
does in H 2 . The relative well depths with angular 
orientation should provide another measure of the 
kind of one-particle resonance involved in low energy 
electron excitation of vibrations. Contrary to H 2 
there will probably be a crossing for coline^r as well 
as C 2v collisions which could, in part, account for the 
higher N 2 cross section [17]. The character of the 
negative-ion should change as a function of angle 
and as a result so should the vibrational-electronic 
coupling. 
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